On the orbit of the short-period exoplanet WASP- 19b 



Coel Hellier\ D.R. Anderson\ A. Collier Cameron^, G.R.M. Miller^, D. Queloz^, B. 
Smalley\ J. Southworth^ & A. H.M.J. Triaud^ 



O 
(N 

D 

(N 
(N 

6 

(N 
> 
m 
a^ 

(N 
m 



ABSTRACT 

WASP-19b has the shortest orbital period of any known exoplanet, orbiting at only 1.2 times 
the Roche tidal radius. By observing the Rossiter-McLaughlin effect we show that WASP-19b's 
orbit is aligned, with A = 4.6 ± 5.2°. Using, in addition, a spectroscopic usin/ and the observed 
rotation period we conclude that the obliquity, ip, is less than 20°. Further, the eccentricity of 
the orbit is less than 0.02. We argue that hot Jupiters with orbital periods as short as that of 
WASP-19b are two orders of magnitude less common than hot Jupiters at the 3-4-d 'pileup'. We 
discuss the evolution of WASP-19b's orbit and argue that most likely it was first moved to near 
twice the Roche limit by third-body interactions, and has since spiralled inwards to its current 
location under tidal decay. This is compatible with a stellar tidal-dissipation quality factor, Q'^, 
of order 10^. 

Subject headings: stars: individual (WASP-19) — planetary systems 



1. Introduction 

Of the known transiting exoplanets WASP-19b 
has the shortest orbital period, at only 0.79 d, and 
thus is important for constraining models of the 
orbital evolution of hot- Jupiter planets. This area 
of theory is currently undergoing rapid develop- 
ment owing to the finding of planets with appar- 
ently short dynamical lifetimes owing to strong 
tidal interactions with their host stars (e.g. Ra- 
sio et al. 1996; Levrard, Winisdoerffer & Chabrier 
2009; Matsumura, Peale & Rasio 2010; Penev 
& Sasselov 2011), of which WASP-18b (Hehier 
etal. 2009) and WASP-19b (Hebb etal. 2010) are 
among the most extreme. 

Further, many of the orbits of hot Jupiters are 
being found to be misaligned with the stellar spin 
axis, with some being retrograde (e.g. Winn et al. 
2009; Narita etal. 2009; Triaud etal. 2010), which 
cannot be accounted for solely by the 'migration' 
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scenario (e.g. Lin et al. 1996) in which hot Jupiters 
form further out and migrate through a protoplan- 
etary disc to short-period orbits. 

Thus a picture is emerging in which hot- Jupiter 
orbits result from a mixture of processes includ- 
ing: disc migration; planet-planet scattering and 
the Kozai mechanism, by which planets can be 
driven into eccentric or misaligned orbits through 
the infiuence of a third perturbing body; and tidal 
dissipation, which circularises and shortens orbits 
(e.g. Fabrycky & Tremaine 2007; Nagasawa, Ida 
& Bessho 2008; Matsumura et al. 2010; Naoz et al. 
2010). 

The relative importance of these mechanisms 
for the ensemble of hot Jupiters can be inves- 
tigated particularly by using the transiting exo- 
planets (e.g. Triaud etal. 2010; Morton & John- 
son 2011), since the Rossiter-McLaughlin effect (a 
distortion of the star's line profiles as the planet 
transits the stellar disc) tells us the angle between 
the planetary orbit and the sky-projection of the 
stellar spin axis. 

Ford & Rasio (2006) pointed out that the or- 
bital period distribution of hot Jupiters cuts off 
near 2 Roche radii, which could be explained if 
planets are scattered into highly elliptical orbits 
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from much further out and then circularise while 
conserving angular momentum (in contrast, disc 
migration alone would predict a smooth distribu- 
tion down to the Roche limit). 

WASP- 19b is one of a small number of plan- 
ets with substantially smaller orbital radii, and 
is the most extreme at only 1.2 times the Roche 
limit. Guillochon, Ramirez-Ruiz & Lin (2010) ar- 
gue that it would be hard to scatter a planet into 
such an orbit, since it would instead by destroyed 
or ejected, and argue that either disc migration 
prior to scattering or tidal decay of the orbit after 
scattering is also necessary. 

We report here an observation of the Rossiter- 
McLaughlin (R-M) effect for WASP-19b, which 
we analyse along with other data to further con- 
strain WASP-19b's orbit, and discuss how this 
planet fits into the above theoretical picture. 

2. Observations 

Wc obtained 25 spectra of WASP-19 through 
a transit on 2010 March 19 using the HARPS 
spectrograph on the ESO 3.6-m at La Silla. The 
data were reduced using the Data Reduction Soft- 
ware by removing the blaze function and cross- 
correlating with a G2 mask (see Triaud et al. 2010 
for further details). In addition 2 HARPS spectra 
were taken two days earlier and 3 more over the 
next two days to better constrain the in-transit 
data. We also report 6 more HARPS spectra and 
3 spectra from the CORALIE spectrograph on the 
Swiss Euler 1.2-m telescope obtained to tie down 
the eccentricity of the orbit (the 39 new radial ve- 
locities are listed in Table 1). In the following 
analysis we also include the 34 CORALIE radial 
velocities previously listed by Hebb ctal. (2010). 

To constrain the transit time, duration and 
ingress and egress lengths, thus improving the 
R-M analysis, wc included a photometric Gunn r- 
band transit obtained using the ESO NTT/EFOSC 
on the night of 2010 February 28, 19 days from the 
HARPS transit. The telescope was dcfocused so 
that each PSF covered approximately 8500 pixels. 
The data were reduced using aperture photometry 
relative to an optimal combination of six compar- 
ison stars following the methods of Southworth 
etal. (2009). We also included 17162 WASP- 
South photometric data points from 2006 -2008 
and the FTS 2;-band lightcurve reported by Hebb 



Table 1: Radial velocity measurements of WASP- 
19. The first three are from CORALIE and the 
remainder from HARPS. 
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etal. (2010). 

The time of secondary occultation (of the 
planet by the star) helps to tie down the eccen- 
tricity, and thus we included in our analysis the 
iJ-band occultation photometry reported by An- 
derson et al. (2010), detrending it with a quadratic 
function of time and sky-background, as in that 
paper, and also the if-band data of Gibson et al. 
(2010), with a linear detrending in time. 

3. Analysis 

To model WASP-19b's orbit we used a Markov 
Chain Monte Carlo (MCMC) approach similar 
to that used by Triaud et al. (2010) and based 
on an R-M model code described more fully 
in Triaud etal. (2009). The datasets listed in 
the previous section were fitted simultaneously 
with a model based on the parameters Tc, P, 
A f , Ti4 , b, Ki Tcff, [Fe/H], ^ecoscj, ^esino;, 
sin /cos A, V v sin / sin A, and, for the occulta- 
tion data, Ai^i.g^m and Ai^2.09^m- Here, Tc is 
the epoch of mid-transit, P is the orbital period, 
AF is the fractional flux-deficit that would be 
observed during transit in the absence of limb- 
darkening, Ti4 is the total transit duration (from 
first to fourth contact), b is the impact param- 
eter of the planet's path across the stellar disc, 
Ki is the stellar reflex velocity semi-amplitude, 
Teff is the stellar effective temperature, [Fe/H] is 
the stellar metallicity, e is the orbital eccentricity, 
w is the argument of periastron, and Ai^i.e/im 
and Ai^2.09/im arc the planet/star flux ratios at 
1.6 //m and 2.09 /im respectively. The resulting 
system parameters are listed in Table 2 and the 
model fits are illustrated in Figure 1. 

As inputs to the MCMC we used Gaussian 
"priors" of T^b = 5500 ±100 K and [Fe/H] = 
0.02 ±0.09 (from Hebb etal. 2010). We also ap- 
plied a prior on vsinl of 5.0 ±0.3 kms~^ (where 
we use vsinl for the projected stellar rotation to 
distinguish from i as the inclination of the planet's 
orbit). This was estimated from the HARPS spec- 
tra by fitting the profiles of several unblended Fe i 
lines. We measured an instrumental broadening of 
0.06 A from telluric lines and we assumed a value 
for macroturbulence of 1.7 ±0.3 kms~^ (Bruntt 
etal. 2010). 

As discussed by Triaud etal. (2011) the choice 
of V sin / prior is critical in modelling the R-M ef- 
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Fig. 1.— (Top) The HARPS (brown) and 
CORALIE (blue) radial velocities of WASP-19 to- 
gether with the fitted model. (Second panel) The 
transit region shown expanded. (Lower panel) The 
NTT transit lightcurve and fitted model. 
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feet when the planet has a low impaet parameter. 
This is because an orbit at right-angles to the spin 
axis and with zero impact parameter is insensitive 
to V sin/. However, for WASP-19b the impact pa- 
rameter (b = 0.66 ± 0.02) is high enough to remove 
the degeneracy, and omitting the prior makes lit- 
tle difference (the R-M angle. A, is 4.6 ± 5.2° with 
the vsinl prior and 6.2 ±6.8° without). 

The stellar Tes, along with the stellar density 
and metallicity, was propagated to the stellar mass 
using the calibration of Enoch etal. (2010). Tes 
was also used to select limb-darkening coefficients 
from the four-parameter models listed by Claret 
(2000; 2004). We also added jitter of 14.1 ms"^ 
to the CORALIE radial velocities and 6.9 ms~^ 
to the HARPS radial velocities in order to obtain 
a spectroscopic xt ^ during the fit and thus 
balance the different datasets in the MCMC. 

The resulting fit (Fig. 1) has a of 176 {u = 

68) for the radial-velocity data (calculated with- 
out added jitter), which likely indicates that stellar 
activity is affecting the radial-velocities, and that 
there are non-uniformities on the face of the star 
during the transit. Note that WASP- 19 has shown 
star spots during transits (Southworth, unpub- 
lished data) , that the star shows a rotational pho- 
tometric modulation at 10.5 d (Hebb etal. 2010), 
and also that the most discrepant HARPS point 
during transit (at phase 1.032) was taken 31 days 
prior to the main HARPS transit. Further, the 
size of the residuals is in line with that expected 
from activity in spotted G stars (Saar & Donahue 
1997). The fit to the NTT transit lightcurve has 
ax^ of 69 {u = 70). 

4. Results and Discussion 

We find that WASP- 19b is in an aligned orbit, 
with A = 4.6 ± 5.2°, where A is the angle between 
the planet's orbit and the sky-projected stellar ro- 
tation axis. Further we find that the eccentricity 
is compatible with zero, with a l-cr upper limit of 
0.009 and a 3-cr upper limit of 0.02. 

Hebb et al. (2010) reported a rotational modu- 
lation in the WASP data at a period of 10.5 ± 0.2 
d. With the stellar radius from Table 2 this would 
translate to an equatorial velocity of 4. 77 ±0.13 
kms~^. This estimate of v is compatible both 
with our spectroscopic wsinJ of 5.0 ±0.3 kms~^ 
and with the v sin / fitted to the R-M of 4.6 ± 0.3 



Table 2: System parameters for WASP-19. 



Parameter (Unit) 



Value 



P(d) 
Tc (HJD) 

Tii id) 

Tl2 = T34 (d) 

AF = Rl/Rl 
b 

z(°) 

Ki (km s-i) 
7 (km s~-^) 
ecosw 
esinw 
e 

c.(°) 

'/'mid —occult 

V sin / (km s"-"") 
A(°) 

(Mo) 

R* (R-o) 

log 5* (cgs) 

P* iPe) 
Mp (Mj,p) 

Rp (Rjup) 
logfifp (cgs) 

PP (p.t) 
a (AU) 
Tp,A=o (K) 



0.7888400 ± 0.0000003 
2455168.96801 ± 0.00009 
0.0655 ± 0.0003 
0.0135 ± 0.0005 
0.0206 ± 0.0002 
0.657 ± 0.015 

79.4 ± 0.4 
0.257 ± 0.003 
20.7873 ± 0.0002 
0.0024 ± 0.0020 
0.000 ± 0.005 
0()46+" ™''^ 
< 0.02 (3cr) 
3 ± 70 
0.5015 ± 0.0012 
4.63 ± 0.26 
4.6 ± 5.2 
< 20 (20-) 
0.97 ± 0.02 
0.99 ± 0.02 
4.432 ± 0.013 

1.168 ± 0.023 
1.386 ± 0.032 
3.143 ± 0.018 
0.438 ± 0.028 
0.01655 ± 0.00013 
2050 ± 40 



kms""'^. That match implies that the stellar spin 
axis is nearly perpendicular to the line of sight, 
which suggests that not only is A low but that the 
obliquity itself, ip, is also low. Since the impact 
factor (6 = 0.66) is relatively high, the R-M value 
will sample higher latitudes and thus might be re- 
duced by differential rotation. We therefore take 
the spectroscopic v sin / as a better indicator of the 
projected equatorial velocity. That results in a 2cr 
limit of / > 65°, which translates (using informa- 
tion from Table 2 and equation 9 from Fabrycky 
& Winn 2009) to a 2a upper limit of V < 20°. 

The known hot Jupiter cxoplanets show a 
'pileup' around periods of 3-4 d (e.g. Szabo & Kiss 
2011). The smaller number at longer periods (~ 10 
d) is to some extent a selection effect, resulting 
from the limited time sampling of ground-based 
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transit surveys. The candidate list for Jupiter- 
sized objects from the better-sampled Kepler data 
shows proportionately more at longer periods, but 
still shows a peak at 3-4 d and a decline beyond 
that (Borucki et al. 2011). The reduced number at 
periods below 2 d is definitely real, being seen in 
both ground- and space-based transit surveys and 
in radial- velocity surveys. For example the fre- 
quency of Jupiter-sized Kepler candidates drops 
by more than an order of magnitude below semi- 
major axis ~ 0.03 AU (period ~ 2 d) (Borucki 
etal. 2011), and note that this is a lower limit to 
the fall-off of planets, since as the planets become 
rarer the transit-mimics would become a larger 
fraction of the candidates. 

Similarly, the WASP-South survey, which has 
found the 3 known Jupiter-sized planets with pe- 
riods < 1 d (WASP-18b, Hellier et al. 2009; WASP- 
19b, Hebb etal. 2010; WASP-43b, Hellier etal. 
2011), shows a much reduced 'hit rate' for plan- 
ets below 2 d (Hellier etal. 2010), even though it 
will be most sensitive to such planets, since they 
produce the most transits and also because the 
transit durations are shorter compared to the data 
lengths. Thus the transiting Jupiters with periods 
below 1 d are rare, and found only because the 
WASP survey covers more than an order of mag- 
nitude more stars than Kepler. 

Note, further, that this means that the actual 
frequency of such planets (not just transiting ones) 
will be lower still, since the range of inclinations 
that produces a transit rises steeply for decreas- 
ing semi-major axis for such close-in planets. For 
example WASP-19b would transit for inclinations 
>74°, but, if it were in the pilcup near 3 4 d, 
the limit would be ^85°, a factor 3 reduction 
in probability. Thus a tentative estimate is that 
planets such as WASP-19b are two orders of mag- 
nitude less common than hot- Jupiters at 3-4 d. 
This means that either Jupiter-sized planets can- 
not easily get into such orbits, or that something, 
such as tidal orbital decay, is then rapidly destroy- 
ing them. 

Ford & Rasio (2006) noted that the lower edge 
of the hot- Jupiter pileup was near twice the Roche 
limit, whereas disc migration would be expected to 
produce a smooth distribution down to the Roche 
limit. They argued that this would arise natu- 
rally if many hot Jupiters arrived from much fur- 
ther out by being scattered into highly eccentric 



orbits that were then circularised while conserv- 
ing angular momentum. The finding of many mis- 
aligned and retrograde planet orbits among those 
in the pileup (Triaud etal. 2010; Winn etal. 2010; 
Narita et al. 2010) also argues that these systems 
(or a large fraction of them) arise not from sim- 
ple migration, but from third-body process such as 
planet-planet scattering or the Kozai mechanism 
(e.g. Fabrycky & Tremaine 2007; Nagasawa et al. 
2008; Triaud etal. 2010). 

However, there are 5 hot Jupiters further in 
than 2 Roche radii (see Matsumura etal. 2010), 
of which WASP- 19b is the most extreme. Us- 
ing the parameters of Table 2, WASP-19b, with 
an orbital period of 0.79 d, has a semi-major 
axis of only 1.21 times the Roche tidal radius 
{an « 2.16i?p(M*/Mp)i/3). Guillochon etal. 
(2010) argue that such orbits are unlikely to result 
from the scattering of planets from outside the ice 
line, followed by circularisation, since the planets 
would instead by destroyed or ejected. Thus they 
argue that these planets must have migrated in- 
wards prior to scattering, or have spiralled inwards 
after scattering as a result of tidal decay. 

Given that the tidal decay timescale for WASP- 
19b in its current orbit is likely to be significantly 
shorter than the system age, we suggest, follow- 
ing Matsumura et al. (2010) and Guillochon et al. 
(2010), that the most likely scenario for WASP- 
19b is: formation beyond the ice line; transfer to 
an orbit near 2aR, at the short-period edge of 
the hot- Jupiter 'pileup', by scattering or by the 
Kozai mechanism, followed by circularisation at 
2 Or,; and then tidal decay of the orbit to the cur- 
rent 1.2 fflR . 

The timescale for tidal decay of the orbit is set 
by the tidal dissipation in the star, denoted by 
the stellar quality factor Q'^, while the eccentric- 
ity damping is also affected by the planetary dissi- 
pation, Q'p, and so could proceed faster, allowing 
circularisation at 2 or before significant decay of 
the orbit (Matsumura etal. 2010). To investigate 
whether this scenario is consistent with the ob- 
served parameters of WASP- 19b we need a value 
of Q'^ small enough such that WASP-19b's orbit 
will decay from '^2aj{ to 1.2 within the life- 
time of WASP-19, estimated at ^ 1 Gyr by Hebb 
etal. (2010), yet large enough to result in a long- 
enough lifetime of WASP- 19b to give a reasonable 
probability of now observing it. 
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Using eqn 5 of Levrard et al. (20090 we find 
that to have decayed inwards from 2 or within 1 
Gyr would require a Q'^ no higher than 10^, while 
a 10-Gyr age would allow up to 10® (though note 
that we also need time for the third-body inter- 
actions leading to the starting point of ^2aR). 
Using the same equation the remaining lifetime 
would be 40 Myr for Q'^ = lO'^, thus giving a^A% 
probability of catching the planet in its current 
state. This probability is in line with the much 
smaller number of Jupiters below 1 d compared to 
at 3-4 d. 

While the range Q'^ = 10^-10® is much higher 
than values taken from binary stars, Penev & Sas- 
selov (2011) argue that, since stars are not spun 
up by planetary-mass companions, the tidal forc- 
ing by planets does not resonate with stellar tides, 
and thus that dissipation is much lower than in 
binary stars, and hence Q'^ can be as high as 10^. 

According to Matsumura et al. (2010) the obliq- 
uity would be damped on a similar timescale to 
the orbital decay, and the eccentricity would be 
damped on either the same or a faster timescale 
(depending on Q'p). Thus this scenario of signif- 
icant orbital decay is consistent with our finding 
of A, ip and e compatible with zero, even if these 
were previously larger during the evolution to and 
at 2 Or . Winn et al. (2010) have argued that plan- 
ets around cool stars with significant convection 
zones (Tcff < 6250 K) will rapidly evolve to aligned 
orbits, and this fits with WASP-19's T^a of 5500 
K. 

Note that it is harder to find a plausible sce- 
nario if we start the infall from larger orbital sep- 
arations. Given the steep dependence of timescale 
on a (r oc (a/i?*)^; Levrard etal. 2009), a small 
enough Q'^, to fit within the age would then pro- 
duce a very small probability of finding it in its 
current orbit. For this reason, it is hard to explain 
WASP-19b's orbit by disc migration and tidal de- 
cay alone: if these processes were efficient enough 
to have moved WASP-19b that far inwards, the 
same processes would likely have destroyed it. 

The other possibility outlined by Guillochon 
etal. (2010) is that WASP-19b first migrated in- 
wards to well within the ice line and was then scat- 
tered to its present semi-major axis. This requires 



^See Matsumura etal. (2010) for corrections to some of the 
equations in this paper. 



a lot of fine tuning to produce a as small as 1.2 or. 
The maximum radius at scattering can be only 
0.01 of the ice line (Guillochon etal. 2010), the 
eccentricity after scattering must not be too high, 
and one has to avoid destruction by tidal orbital 
decay; yet if there hasn't been significant tidal de- 
cay then one could expect relics of the scattering 
in the form of eccentricity or misalignment. Thus 
our findings of an aligned, circular orbit make this 
already unlikely scenario even less likely. 

Thus, the current orbit of WASP-19b, the tight- 
est of the known exoplanets, is most plausibly ex- 
plained if the planet first moved to near ^ 2aR by 
scattering or by the Kozai mechanism, as has been 
suggested for many of the hot Jupiters in the 3-4- 
d pileup, and circularised there; and then, perhaps 
through starting at the lower edge of the pileup, 
spiralled inward to its current location by orbital 
decay. 

Based on observations made with the ESO 3.6- 
ni/HARPS, program 084-C-0185, and the ESO 
NTT, and the Euler 1.2-m at La Silla Observa- 
tory. 
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